Topological Origin of Zero-Energy Edge States in Particle-Hole Symmetric Systems 
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A criterion to determine the existence of zero-energy edge states is discussed for a class of particle- 
hole symmetric Hamiltonians. A "loop" in a parameter space is assigned for each one-dimensional 
bulk Hamiltonian, and its topological properties, combined with the chiral symmetry, play an essen- 
tial role. It provides a unified framework to discuss zero-energy edge modes for several systems such 
as fully gapped superconductors, two-dimensional d-wave superconductors, and graphite ribbons. A 
variants of the Peierls instability caused by the presence of edges is also discussed. 



Depending on several parameters such as hopping in- 
tegrals or chemical potentials, and also on underlying 
crystalline lattices, a large variety of electronic structures 
are realized in condensed matter physics. Electron cor- 
relations also give rise to a plenty of quantum phases, 
forming non-trivial quasi-particle band structures. An 
interesting consequence of a rich band structure is the 
existence of edge states that may appear when bound- 
aries are present. In the quantum Hall effect (QHE), this 
issue was discussed in terms of the origin of the quanti- 
zation of a Hall conductance. ||, ^, || ^ Recently, 
the ideas developed in the QHE have also been extended 
for other gapped many-body systems, and become essen- 
tial to describing topological nature of several quantum 
phases. |, 0, 1 H 

Apart from these examples for gapped systems, edge 
states in gapless systems have attracted much attention 
recently. Example of these are d-wave superconductor 
(SC) with edges 14|, or graphite ribbons where 
the existence of edge states strongly depends on the shape 
of edges. For d-wave SC with edges, the zero bias conduc- 
tance peak (ZBCP) due to zero-energy ed ge states was 
observed via a tunneling spectroscopy, fl^, [T^ 

The issue addressed in this Letter is how to infer the ex- 
istence of zero-energy eigen states localized on the bound- 
aries in terms of properties of the bulk, and the sym- 
metry. We first consider one-dimensional (ID) systems 
with a particle-hole symmetry, and then apply the re- 
sults to systems in higher dimensions. Especially, we will 
demonstrate applications to fully gapped SC in conjunc- 
tion with the Chern number, 2D d-wave SC, and graphite 
ribbons. In addition to these examples, the present work 
is also applicable to zero-modes in the ID molecule poly- 
acetylene [Q, and quantum spin systems. 

We start with the following single-particle Hamiltonian 
on a ID lattice: 



A' 
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where ts;^^' , A^^^i , A'^ ,^, e C, and cj, = (c^pC^;;) denotes 
electron creation/annihilation operators at site x. The 
total number of the lattice sites is and x — 1, ■ ■ ■ , N^- 
This Hamiltonian includes the Bogoliubov-de Gennes 



(BdG) Hamiltonian both for singlet and ( some of ) 
triplet SC. 

In the following, we consider two types of Hamilto- 
nians: the bulk Hamiltonian and the edge Hamiltonian. 
As for the bulk Hamiltonian, assuming that the system is 
translationally invariant, hx,x' = h{x — x'), and adopting 
periodic boundary condition (PBC), we can perform the 
Fourier transformation to obtain 7^'^"''' = cj^hkCf, = 
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Cj,, where = 1/a 



Afe 

-a 

= is the crystal momentum, and ^fc € 
Since {(TYhkCY)* = —hk {o'x,y,z represent 
the Pauli matrices), eigenvalues E and —E always appear 
in pair for each fc, which we call the particle-hole sym- 
metry. Let us introduce a convenient parametrization for 
^buik k-space 

hk — R{k) ■ cr, 

where R{k) = {X,Y,Z) := (ReAfe, -ImA^, ^fc) e R^. in 
this parameterization, the energy eigenvalues are given 
by E{k) = ±\R{k)\. The origin O eM? corresponds to 
the gap-closing point. For a given k £ , there exists a 
one to one correspondence between a point in a 3D space 
R(k) and hk, and hence we can identify a loop £ : k € 
Rik) e R3 for each ID Hamiltonian T^^"'"^; for a 
given parametrized loop, we can always reconstruct Ti'^""' 
by inverse Fourier transformation. We write 7i''""'[£] 
for the Hamiltonian which corresponds to £ hereafter. 
An edge Hamiltonian 7i°'^so jg generated by truncating 



a bulk Hamiltonian Ti. 



bulkr 



in a certain way. We refer 



an edge Hamiltonian as Jf'^'^SB gj ^ -^^here e represents a 
prescription for creating edges. For example, a natural 
way of truncation is to prohibit all the matrix elements 
across N^, i.e., set hx,x' = ii & [x, x'], which we call 
Cc. Generally, e can represent an impurity potential at 
an edge, coexistence of different order parameters near 
boundaries in superconducting systems, etc. Then, we 
ask if W'^so[£ g] supports zero-energy states localized at 
either end of the sample for given £ and e. Our strategy 
to answer this question is to consider a continuous defor- 
mation of a Hamiltonian from a reference Hamiltonian 
with exact zero-energy edge states, in conjunction with 
a symmetry. 
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In the following, let us focus on a loop on a 2D plane 
that contains the origin O in i?-space. We crown such 
loops with a superscript * as a reminder, thereby referred 
to as £* . As a prescription for creating edges, we adopt 
Cc for a while. Let | £,E,p) denote an edge states of 
■^edgej^^ ec] with energy E, localizing at p = L{R) where 
L{R) represents the left (right) edge. We assume a state 
which appears within the bulk energy gap is localized 
at either end of the sample for an infinite system. A 
state localized at both ends also may appears, which is 
a superposition made from two independent edge states 
localized at the left and right. We will show 

(A) if H'='^so[£*^gj has an edge state at 
non-zero energy | £*, i? 0,p), it also has 
I i*,—E,p) which locahzes at the same edge, 
with the opposite energy. 

First, note that we can restrict ourselves to loops on 
the A"F-plane, since an arbitrary 2D plane can be rotated 
to the XY-plane by a unitary transformation: a global 
50(3) rotation in il-space, which amounts to a SU{2) 
transformation on for each site. To prove the state- 
ment, it is essential that the particle-hole symmetry is 
promoted to the chiral symmetry for ']-C'^'^s° [£* , ■ Since 
all the hopping tx/x is zero for loops on the XY-plane, 
the Hamiltonian can be expressed as Ti'^'^so j£* ^ g^j — 



(ct , 



)H 



H 
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anticommutes with YRY = —H, which we call the 
chiral symmetry. Consequently, if Ji'^'^sc j^* ^ g^j j^^s an 
edge mode | ip) = \ £* , E ^ 0,p), it also has an edge mode 
with energy —E, Fj V') = | £*,—E,q). Moreover, since F 
is a purely local operator which only changes the phase 
of c-f, i.e., it does not "mix" the coordinate in the real 
space, I ip) and F| ip) should be localized at the same 
edge, q = p. Notice the above discussion is not applica- 
ble for E ^ 0, since both | ip) and F| ip) have the same 
energy, and hence can be the equivalent state. 

Next, we further assume that £* is continuously de- 
formed into a unite circle £c centered at O, such that the 
loop is always on the 2D plane, and does not cross O 
during the deformation. (Fig. [l]). For a loop £* with this 
property, we write as £* £c henceforth. We can prove 
that: 

(B) 'W^^^[£* ~ £c,ec] has at least a pair of 
edge states at zero energy. 

To see this, we focus on R{k) = (cos fc, — sin A:, 0) 
and the corresponding Hamiltonian Ti'^^scj^^^ g^j _ 
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Cx+i + h.c. Since c^, cj^, Cj^^^ cj^^^ do 



not appear in W'^scj^^^ g^j^ there are two exact zero- 
energy levels, which localize at a; = 1 and x — Nx, i.e., 
n'"^s<^[£c, ej has two edge states | 4, 0, L) and | £c, 0, R). 
pO| By assumption, we can deform £c into £* continu- 
ously. During the deformation, | £c,0,L) and | £c,0,R) 
do not go away from zero energy, since we can apply (A) , 




FIG. 1: (a) Continuously deforming £c into a loop £* ~ £c. 
During the deformation, the loop is kept on the 2D plane, 
without crossing O. (b) A possible energy spectrum during 
the deformation. A thick/broken line represents a edge mode 
localized at R/L. 



and the bulk energy gap does not collapse. Although 
other edge states | E,p) and | —E,p) may appear in pair 
from the bulk energy bands, since the number of edge 
modes localized at L/R is always odd, there must exist 
at least a pair of zero-energy states. 

Although we have concerned ourselves with a certain 
type of edges Cc, let us next consider to adiabatically 
modify Cc- As far as the modified prescription does not 
break the chiral symmetry, the perturbed Hamiltonian 
also supports exact zero-energy edge states, since per- 
turbations at the edges do not collapsed the energy gap. 
Thus, we have showed 

(C) for a prescription e* that respects the chi- 
ral symmetry, W^scj^* ^ £c,e*] possesses at 
least a pair of zero-energy states. 

In summary, there are three conditions for H'^'^^'^[£,e] 
to support zero-energy edge states: (A) £ is on a 2D plane 
that contains O, {£*), (B) £ is continuously deformed to 
£c without crossing O, {£ ^ £c), and (C) e respects the 
chiral symmetry, (e*). 

We have established our main results, and a few com- 
ments are in order. First, notice that the edge states 
discussed here are not at exact zero energy for a finite 
system, though | £c,0,L) and | 4,0, i?) are exact zero- 
energy state. This is allowed since an assumption for 
the statement (A) does not hold for a finite system size. 
In this case, a state localized at both ends cannot be 
decomposed into two-independent edge states, which we 
can regard as a hybridized state made from the two edge 
modes at the left and right. In Nx — + oo, this state be- 
comes degenerate with another hybridized state. 

Second, consider a unit circle that encloses O n 
times (n:odd). 'H°'^^°[£^,ec] can be diagonalized in the 
same way as Tl°'^^°[£c, ed , resulting in 2n exact zero- 
energy states. Then, by the same discussion, a class of 
Hamiltonians H^'^^^[£* ~ ^",e*] have at least one pair of 
edge states at E = 0. 

Finally, the present discussion is consistent with 
the King-Smith- Vanderbilt (KSV) formula that relates 
macroscopic polarization to the Zak's geometric phase 7. 
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FIG. 2: Loops in Ji-space and the energy spectrum of d^2_y2- 
wave SC with (a) (110) and (b)(100) surfaces. Dotted squares 
show a choice of unit cell in Fourier transforming along the 
edges. The calculation is for Nx = 50 for (110) surfaces, and 
iV:^, = 30 for (100) surfaces. 



|2|, |4| g| However, the KSV formula does not tell 
us if edge states are at zero energy. 

We go on to applications of the present results. We 
adopt Cc as a prescription for creating edges unless oth- 
erwise stated. First, we discuss fully gapped systems and 
edge states. Especially, we comment on a topological as- 
pect of 2D SC with a full gap, whose examples include 
d + id SC and the chiral p-wave SC. |, 0, |, |, |o[ |l|] 
For these SC, we can define an integer called the Chern 
number, non-zero value of which implies the existence of 
edge states connecting the upper and the lower bands as 
known in the QHE Q . The present results are consistent 
with this discussion. For 2D systems with edges, we first 
Fourier transform along a direction parallel to the edge, 
to get a family of ID Hamiltonians parametrized by the 
wave number along the edge. Then, we can apply the 
present discussions for each ID Hamiltonian. Since the 
non-zero Chern number implies there exists a loop which 
is on a plane and encloses O [0, both the topologi- 
cal argument and the present results lead to existence of 
zero-energy edge modes. For fully gapped systems, edge 
modes are expected to be stable even in the presence of 
electron-electron interaction as far as the bulk energy gap 
is not collapsed. 

Although the topological argument is only applicable 
for fully gapped systems, our results here are not re- 
stricted to gapped cases, and can be applicable also for 
gapless cases in arbitrary dimensions. Here, as an ap- 
plication, we consider surface states for dj.2_j,2-wave SC. 
In Ref. ||l^, a semi-classical approach was employed to 
show the sign change of the pair potential at a (110) sur- 
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FIG. 3: Loops in it-space and the energy spectrum of a 
graphite ribbon with (a) zigzag, (b) bearded, and (C) arm- 
chair edges. The ovals indicate how to form a spinor c, for 
each edge, and dotted squares show a choice of unit cell 
in Fourier transforming along the edges. The loops corre- 
sponding to a one-parameter family of Hamiltonians are (a) 
Rkyik^i) = (cos(fcj; — fc^/) + 1 + cosfcy, — sin(fcy — k,.,) + 



sinfcj;,0), (b) (cosfca;' + cos{ky 



k-i- 



+ 1, sin k^i — sin(fcy — 
fc^/),0), (c) (cos(fcj./ -I- ky) + cosfcj;/ + 1, — sin(fc^/ + ky) + 
sinfcjj/jO). Here we have taken all the hopping integral equal 
to unity, and ky is a wave number along the edges. The cal- 
culation is for N^i — 30 for zigzag and bearded edges, and 
N^i = 29 for armchair edges. 



face gives rise to existence of edge states, which can be 
used as a phase sensitive probe to detect pairing symme- 
tries. It was also pointed out the Andreev equation for 
the present system is closely related to Witten's super- 
symmetric quantum mechanics. [|3[ |l^ Here, we discuss 
this issue with a lattice regularization. 

Consider 2D dx2_ 

^PBC 



[clh. 



, hy — 



wave 
■ h.c 
t -A 
-A -t 



SC H^""^ 
f ct/iQC^] , where 
' /i 



and /iQ = 



-^l 



t A 
A -t 

(We set i = A = 1, /i = as an example.) We terminate 
this system, and consider (110) surfaces first. Fourier 
transforming along the y' direction in Fig. |^ (a), 
we obtain a family of ID Hamiltonians parametrized 
by kyi . The corresponding loops are ,{kx) = 



4 



(2 cos(fca; — kyi) — 2 cos kx, 0, 2 003(^3; — fcj,/) + 2 cos kx)- 
Foragiven ky, , {l+cosky,){X/2f + {1- cos ky,){Z/ 2)^ = 
2 sin^ fcj,' is satisfied, which is an ellipsis on the XZ- 
plane enclosing O. Thus, from the above discussion, 
the present system supports zero-energy surface states 
for all kyi except at the gap-closing points kyi — ztn, 
where the loop collapses into a line segment. 

On the other hand, for (100) surfaces, we obtain 
Rky{kx) = (2(cosfcx — cos fcj^), 0,2 (cos fc^c-f cos fcy)), which 
is a line segment on the XZ-plane for all ky. Zero-energy 
edge states are not expected to exist for this case. We 
have verified numerically this prediction in Fig. |^(b). 

Let us comment on an interplay between zero-energy 
edge states and interactions for the present case. If we 
treat the problem self-consistently, coexistence of is- or 
idxy-'Vf&'ve order parameter with d2,2_j^2-wave near the 
surface is possible for the (110) surface, locally breaking 
the time-reversal symmetry, pq] . This can be interpreted 
based on the present discussions as follows. Since edge 
states with different ky> are all degenerate at £' = 0, 
they are expected to cause a Peierls-like instability. In 
presence of interactions, parameters in a single particle 
Hamiltonian t, A, A' near the edges might be effectively 
modified in order to lift the degeneracy, and thereby 
lower the ground state energy. However, since these zero- 
energy edge states are stable to perturbations which re- 
spect the chiral symmetry (statement (C)), such modifi- 
cations should be accompanied with the breaking of the 
chiral symmetry near the boundaries. The emergence of 
is or idxy components near the boundary indeed breaks 
the chiral symmetry to lift the degeneracy of edge modes, 
while a purely real order parameter cannot do it. 

We turn to edge states in graphite ribbons. There 



are several types of edges for a graphite ribbon such as 
zigzag, bearded, and armchair edge.|p^ Defining C| = c, 

and c| — Co, where c,/o is an electron annihilation oper- 
ator on a sublattice •/o, we can apply our formalism 
to graphite ribbons. Notice that we have several op- 
tions for choosing c,/o to form a spinor c, since they 
live on different sites. When we truncate the system, 
these choices lead to different shapes of edges (Fig. ||). 
Taking an appropriate pair for each type of edge as indi- 
cated in Fig. H, we can discuss in parallel to the above 
SC example. The existence of zero-energy edge states is 
predicted for the zigzag and the bearded case, while we 
do not expect zero-energy edge states for an armchair 
edge, which is confirmed by a numerical calculation (see 
Fig.y). These zero-energy edge modes are continuously 
connected to the gapless bulk spectrum, forming a flat 
band and a sharp peak in density of states at the fermi 
energy. This might trigger an instability in presence of 
electron-electron or electron-phonon interactions, which 
leads to, for example, a magnetic polarization near the 
boundaries, jl^ 

To conclude, we have established a criterion to deter- 
mine the existence of zero-energy edge modes in terms of 
bulk properties and the chiral symmetry. Our strategy 
is to make use of the chiral symmetry, and a continu- 
ous deformation of a reference Hamiltonian with exact 
zero-energy edge states. The present discussions are ap- 
plicable for both gapped and gapless systems in arbitrary 
dimensions. 
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